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SUMMARY 


Uf\  - 


^ The  orbit  of  Proton  4,  1968-103A,  -has-beejj- redetermined,  in  greater  detail 
and  with  better  accuracy,  -in-order5  to  clarify  previously  puzzling  features  in  the 

'rr-c  ' 

variation  of  orbital  inclination.  Orbital  parameters  have-been^determined  at 
25  epochs  between  December  1968  and  July  1969,  using %about  1600  optical  and  radar 
observations  with^the  RAE  orbit  refinement  program  PR0P6. 

During  January  1969  the  orbit  passed  through  31:2  resonancejj^when  the 

ground  track  over  the  Earth  repeats  every  two  days  after  31  revolutions  of  the 

satellite.  A simultaneous  least-squares  fitting  of  theoretical  curves  to  the 

values  of  inclination  and  eccentricity  bet  'een  14  December  1968  and  6 March  1969 

-has^yielded  values  for  two  pairs  of  lumped  31 st-order  geopotential  coefficients, 

A >■  V'vf?  # t t 

appropriate  to  an  inclination  of  51.5  ^‘'This  -i-sx  the  first  specific  evaluation  of 
31 st-order  coefficients. 

The  15  values  of  inclination  after  the  resonance,  from  March  to  near  decay 

in  July  1969,  have-bl^n^used  to  determine  mean,  morning  and  afternoon-evening 

values  for  the  rotation  rate  of  the  atmosphere  at  a height  near  260km;  the  values 

CU 

of  rotation  rate,  namely  1.1,  0.9  and  1.3rev/day  respectively,  confirm  '•the''  trends 
already  established  from  analysis  of  other  satellite  orbits. 
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1 INTRODUCTION 

Proton  4,  1968-103A,  was  launched  on  16  November  1968  into  an  orbit  of 
inclination  51.5°;  initially  the  perigee  height  was  near  250km  and  apogee  near 
480km.  The  satellite  was  extremely  massive,  weighing  17000kg,  and  was  in  the 
form  of  a squat  cylinder  3m  long  and  4m  in  diameter;  it  remained  in  orbit  for 
250  days  and  decayed  on  24  July  1969. 

The  orbit  seemed  promising  for  use  in  geophysical  studies,  and  was  deter- 
mined* in  1970  at  20  epochs  during  its  life,  using  US  Navy,  optical  and  radar 
observations.  Analysis  of  the  decrease  in  inclination  yielded  an  acceptable 
value  of  the  atmospheric  rotation  rate,  but  there  were  some  peculiarities  in  the 
variation  of  inclination,  particularly  a 'dip'  between  mid  January  and  early 
March  1969.  It  was  suggested  in  1970  that  this  might  be  due  to  neglected  effects 
of  tesseral  harmonics  in  the  geopotential.  The  orbit  has  now  been  redetermined 

at  25  epochs,  from  about  1600  US  Navy,  optical  and  radar  observations,  using  the 
2 3 

improved  version  ’ of  the  RAE  computer  program  for  orbit  refinement,  PR0P6,  and 
assigning  more  realistic  a priori  accuracies  to  the  US  Navy  observations. 

The  new  orbit  determination  shows  that  the  'dip'  in  inclination  was  the 
result  of  a strong  resonance  effect  associated  with  tesseral  harmonics  in  the 
geopotential  of  order  31,  at  the  time  when  the  ground  track  of  the  satellite 
over  the  Earth  was  repeating  every  two  days,  after  31  revolutions.  Since 
Proton  4 was  so  massive,  it  was  less  affected  by  drag  than  most  satellites  at 
similar  heights,  and  passed  through  the  resonance  slowly  enough  to  allow  values 
of  lumped  31st-order  geopotential  coefficients  to  be  determined.  Values  of 
atmospheric  rotation  rate  are  also  obtained  for  the  last  four  months  of  the 
satellite's  life. 

2 ORBIT  DETERMINATION 
2. 1 Observations 

More  than  1600  observations  were  available,  made  between  December  1968  and 
July  1969.  In  the  previous  orbit  determination,  all  the  US  Navy  observations 
were  used,  but  in  this  redetermination  of  the  orbit  165  US  Navy  observations 
having  elevations  less  than  20°  were  omitted,  because  they  are  known  to  be  of 
poorer  accuracy  than  observations  at  higher  elevations.  Many  of  the  US  Navy 
observations  are  expressed  relative  to  the  Earth's  centre,  and  the  a priori 
angular  accuracy  assigned  to  these  observations  has  been  changed  from  0.02°  to 
0.01°.  Six  observations  were  used  from  the  Hewitt  camera  at  Malvern,  seven  from 
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the  200mm  camera  at  Meudon,  and  24  from  the  kinetheodolite  at  the  South  African 
Astronomical  Observatory.  There  were  about  700  US  Navy  observations,  about  180 
visual  observations  and  about  140  observations  from  the  Malvern  radar.  In  the 
course  of  the  orbit  determinations  about  400  observations  were  rejected. 

2.2  Observational  accuracy 

Table  1 summarizes  the  residuals  obtained  for  stations  contributing  at 
least  seven  observations  to  the  orbit  determination  and  for  the  Malvern  Hewitt 
camera  (with  six  observations).  The  residuals  were  evaluated  using  the  computer 
program  ORES  , and  the  results  have  been  sent  to  the  individua  observers.  The 
most  accurate  observations  are  those  from  the  Malvern  camera  (rms  residual 
4 seconds  of  arc).  The  next  most  accurate  are  the  ceven  from  the  200mm  camera 
at  Meudon  (6  seconds  of  arc).  The  24  observations  from  the  kinetheodolite  at 
the  South  African  Astronomical  Observatory  are  accurate  to  about  40  seconds  of 
arc.  The  US  Navy  observations  (totalling  about  700)  have  topocentric  accuracies 
between  3.4  and  4.5  minutes  of  arc  (station  29  is  geocentric  and  the  residuals 
need  to  be  multiplied  by  a factor  of  about  5 to  obtain  an  equivalent  topocentric 
accuracy,  which  is  3.2  minutes  of  arc).  The  199  Malvern  radar  observations 
(which  include  about  60  removed  from  the  orbit  determinations  to  avoid  bias  on 
occasions  when  radar  observations  were  too  numerous)  have  rms  residuals  of  about 

3.2  minutes  of  arc.  The  180  or  so  visual  observations,  obtained  from  either  the 
Appleton  Laboratory  at  Slough  or  the  (now  disbanded)  Moonwatch  Division  of  the 
Smithsonian  Astrophysical  Observatory,  have  residuals  varying  between  4 and  14 
minutes  of  arc. 

These  residuals  cannot  be  compared  directly  with  those  determined  in  1970, 
since  different  observations  have  sometimes  been  used.  The  US  Navy  observations 
now  appear  more  accurate  than  before,  but  this  is  because  of  the  omission  of 
low-elevation  observations  and  the  reduction  from  0.02°  to  0.01°  in  the  pre- 
assigned accuracy  of  observations  expressed  relative  to  the  Earth’s  centre 
(with  the  result  that  moie  were  rejected).  On  average,  the  visual  observations 
show  little  change  in  their  residuals.  The  abc  residuals  (the  arithmetic  mean 
of  the  best  70  to  80  per  cent  of  the  visual  observations)  are  rather  higher, 
but  differences  are  not  significant  since  the  choice  of  observations  to  be 
eliminated  is  not  an  exact  procedure. 


Table  1 

Residuals  for  selected  observing  stations 


Number 

of 

observations 


rms  residuals 


Minutes  of  arc 


abc  residuals  i 
minutes  of  arc! 


US  Navy 
US  Navy 
US  Navy 
US  Navy 
US  Navy 
29  US  Navy 
719  Rodewisch 
2265  Farnham 

2303  Malvern  (Hewitt  camera) 

2304  Malvern  (Radar) 

2403  Stevenage  3 
2409  Ribbleton* 

2512  Mountcastle** 

2525  Crawley  Ridge 
2531  Ditton  Park  East+ 

2539  Dymchurch 

2573  Genoa  1 

2577  Cape  kinetheodolite 
8030  Meudon  camera 


* Includes  Moonwatch  station  667 
**  Includes  Moonwatch  station  656 
+ Includes  Moonwatch  station  661 

2.3  Orbital  accuracy  * 


The  25  sets  of  orbital  elements  as  determined  by  PR0P6  are  given  in 
Table  2.  The  mean  anomaly  is  represented  in  PROP  by  a polynomial  in  time 
from  epoch, 


M = + M,t  + M„t2  + M„t3  + M,t4  + M„t5 


6 


143 


The  drag  proved  to  be  rather  variable,  and  9 of  thr  25  orbits  required  all  six 
terms,  as  shown  by  values  up  to  Mj.  in  Table  2.  With  the  other  orbits,  the 
observations  were  better  fitted  using  less  coefficients. 

Orbit  1 played  no  part  in  the  final  analysis  because  it  proved  to  be  too 

distant  from  the  31:2  resonance.  Orbits  2-10  were  used  in  the  analysis  of  31:2 

resonance  effects.  The  sd  in  inclination  for  these  9 orbits  varied  from  0.0002° 

(for  orbit  2 with  Hewitt  camera  observations)  to  0.0016°,  the  rms  being  0.0011°. 

“6 

The  sd  in  eccentricity  varied  from  3 x 10  (equivalent  to  about  20m  in  perigee 

“6  “"6 

or  apogee)  to  21  x 10  , with  an  rms  of  14  * 10  . The  corresponding  rms  values 

for  the  first  determination'  made  in  1970,  were  0.0014°  in  inclination  and 

_6 

17  x 10  in  eccentricity  - so  the  accuracy  has  been  improved  by  about  20%.  The 
rms  value  for  e , the  measure  of  fit,  is  0.75  here,  compared  with  0.51  in  1970. 
This  indicates  that  the  pre-assigned  errors  for  the  US  Navy  observations  were 
too  large  previously  and  are  nearer  the  mark  here. 

Orbits  1 1-25  were  used  to  determine  the  atmospheric  rotation  rate:  their 

sd  in  inclination  varied  over  the  range  0.0001°  to  0.0018°,  the  rms  value  being 

o -*6  “6 

0.0011  . The  sd  in  eccentricity  varied  from  10  x 10  to  38  x 10  , with  an  rms 

__6 

of  22  x io  . The  corresponding  1970  mean  standard  deviations  in  i and  e 
were  0.0011°  (no  change)  and  29  x 10  ^ (30%  worse). 

For  all  25  orbits,  the  sd  in  Q varied  between  0.001°  and  0.004°. 

and  io  have  similar  standard  deviations,  varying  between  0.02  and  0.17°  for 

orbits  1-19,  but  increased,  for  orbits  20-25,  to  between  0.18°  and  0.50°.  The 

2 

decay  rate  had  an  rms  sd  of  0.0008deg/day  for  the  first  24  orbits  (an 

accuracy  cf  about  0.2%);  on  the  final  (25th)  orbit,  the  sd  increased  to 
2 

0.0044deg/day  (but  is  accurate  to  0.1%). 

Some  difficulty  was  encountered  in  determining  a few  of  the  orbits.  For 
orbits  12  and  17,  for  example,  magnetic  storms  occurred  near  epoch,  and  observa- 
tions on  the  magnetic-storm  days  had  to  be  removed  before  satisfactory  orbits 
were  obtained.  Orbit  20  was  particularly  difficult  to  determine,  with  both 
Hewitt  camera  and  Meudon  camera  observations  present.  Although  the  Hewitt 
camera  observations  appeared  to  fit  satisfactorily,  they  had  to  be  removed 
before  achieving  an  orbit  which  fitted  in  with  the  other  orbits.  Orbit  24  had 
observations  from  only  two  stations,  mostly  from  a US  Navy  station.  Although  a 
good  orbit  appears  to  have  been  obtained,  it  may  suffer  from  bias  and  should  be 
treated  with  caution. 
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3 ANALYSIS  OF  THE  31:2  RESONANCE 

3. 1 Theoretical  equations 

Near  31:2  resonance,  the  appropriate  resonance  angle  $ is  given  by 

4-  = 2 (to  + M)  + 31  (n  - v)  , (2) 

where  v is  the  sidereal  angle.  Fig.l  shows  the  variation  of  3>  for  Proton  4 
and  indicates  that  exact  31:2  resonance  occurred  on  18  January  1969  (MJD  40239). 
The  resonance  affected  the  orbit  mainly  between  14  December  1968  and  25  February 
1969,  when  <P  was  within  lj  cycles  of  its  value  at  resonance  (318°). 

The  Earth's  gravitational  potential  l at  an  exterior  point  (r,0,A)  may 
be  written"*  in  normalized  form  as 


U = 


U + - 


V V 


LL\r  r 


m 


(cos  0)  I C.  cos  mA  + S„  sin  mA^N 
1 £m  £m  1 


£m 


(3) 


£=2  m=l 


where  U is  the  longitude- averaged  potential,  r the  distance  from  the  Earth's 

centre,  9 the  co-latitude,  A the  longitude  (positive  to  the  east),  p the 

3 2 

gravitational  constant  for  the  Earth  (398601km  /s  ),  and  R the  Earth's  equato- 
rial radius  (6378.1km).  P™  (cos  9)  is  the  associated  Legendre  function  of 

order  m and  degree  £ , and  C.  and  S„  are  the  normalized  tesseral 

£m  £m 

harmonic  coefficients,  of  which  those  of  order  m = 31  particularly  concern  us 
here.  The  normalizing  factor  N is  given  by"*: 


m2  _ 2(2£  + 1 ) (£  - m) 1 

&m  (£  + m) ! 


(4) 


Near  the  31:2  resonance  the  theoretical  variation  of  inclination  may  be 
6,7,8 

written  * ’ 


Ui. 

dt 


sm  1 


[7j  [I  (31  ' 2 C0S  i}  F32,31,15[S31  sin*+C31  C0S  *j 


fj.l  _1.1  1 

+ 1 7e (3 1 - cos  i)  F,^  3]  sin  ($  - <o)  - S3J  cos  (4>-oj)> 


L-1,3  _-l,3  ] 

+ 1 3e (31  — 3 cos  i)  F3J  3]  sin  ($  + to)  - S3)  cos  (<t>  + <o)> 


, _ . I q ! cos  , . . 

+ terms  m e1^'  . (y<J>  - qto) 

s in 


> 


(5) 


In  equation  (5),  y and  q are  integers,  with  y taking  the  values  1,2,3,... 
and  q the  values  0,  ±1,  ±2  ...  . Only  the  terms  with  y = 1 and  q = 0 and 
±1  are  given  explicitly,  since  it  is  believed  that  the  others  will  be  small: 
the  terms  with  |q|  = 2 have  e as  a multiplying  factor  (where  e — 0.015 
for  Proton  4),  while  the  terms  with  y = 2 are  associated  with  harmonics  of 
order  62,  which  should  be  much  smaller  than  those  of  order  31. 

_0,2  J,i 

The  quantities  S3J  , C3J  , etc.  in  equation  (5)  are  what  might  be  called 

'natural  lumped  coefficients'  of  order  31-  In  defining  them  it  is  helpful  to 

7 

consider  the  general  B:ci  resonance  and,  for  given  y and  q , to  introduce 
integer  suffixes  m,  k and  p such  that,  for  suitable  values  of  £ (Ref. 6): 
m = yB;  k = ya  - q ; and  2p  = £ - k . Then  the  lumped  coefficients  may  be 
written^ 


_q>k  V'  _ q,k 

c = ) Q„C„  , S 

m £ £m  m 

£ £ 

where  £ increases  in  steps  of  two  from  a minimum  value  (normally  equal 

to  m or  m + 1)  , and  the  are  constant  factors  with  = 1 . 

We  are  concerned  with  the  31:2  resonance  and  are  considering  only  the 
terms  with  y = 1 ; so  we  have  m = yB  = 31  , and  k = yct-q  = 2-  q . Thus, 
for  the  three  terms  included  in  equation  (5),  with  q = 0,  1 and  -1,  the 
affixes  [q,k]  are  [b,2],  [l , l]  and  [-1,3]  respectively,  as  shown  in  equation  (5) 
The  relevant  values  of  £ , for  use  in  equation  (6),  must  be  such  that  £ > m 

mm  1 , 1 

and  (£  - k)  is  even:  with  for  example,  £ %■  31  and  (£  - 1)  must  be 

even;  hence  £ = 31,  33,  35, 


■l 


Q£°£m 


(6) 


• • • 
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The  quantities  F appearing  in  equation  (5)  are  the  functions  of  inclina- 
tion defined  by  Allan7.  The  factors  in  equations  (6)  depend  on  the 

anc*  °n  > and  ^or  31:2  resonance  the  lumped  coefficients  may  be 

written  explicitly  in  terms  of  the  geopotential  harmonic  coefficients  C as 
follows3: 


_0,2 

C31 


= C 


_ F34 ,31,16  /Rfp  , F36, 31,17  /R  f- 

32,31  F32,31 ,15  a ' 34,31  F32,3I ,15  3 ' 36, 


31 


(7) 


'31 


_-l,3 

C3. 


= C 


31,31 


= C 


31,31 


18F 

33,31 , 16 

fife 

19F 

35,31,17 

fi 

l7F31,3l,15 

\a ) 33,31 

! 7F31 ,31,15 

Va 

14F 

33,31,15 

fife 

15F 

35,31,16 

f- 

13^31 ,31,14 

\a  J 33,31 

13F3l ,31 , 14 

U 

,31 


,31 


(8) 


(91 


and  similarly  for  S , on  replacing  C by  S throughout.  The  numerical  values 
of  the  three  most  important  F functions  used  here  are: 


F32,3I, 15 


0.23668  sin29i  (16  cos  i - I)  (I  ♦ cos  i)2  • 4.6336  * 10  3 for  Proton  4 (i  » SI. 53°) 


r3,  3|  ,5  - 0.49905  sin30i  (I  + cos  -i) 


'31,31,14 


0.44035  sin28i  (I  ♦ cos  i)3 


0.52694  * !0~3  for  Proton  4 


1.9949  « 10  3 for  Proton  4 


(10) 


The  theoretical  equation  for  the  variation  of  eccentricity  near  31:2 

-q,k  -q,k 

resonance,  involves  the  same  lumped  coefficients  C and  S as  the 

6 7 m m 

equation  for  di/dt  , and  may  be  written0 ’ 


de 

dt 


R 

0 U 


31 


R - 

2 

_0,2 

-IF32 

,31, 

15  e 

f3. 

sin 

4> 

+ C^j  cos  $> 

I 

[J,1 

1,1 

17F31 

,31, 

15  ] 

[C3I 

sin 

(0 

- U)  - S31 

1 

L"1, 

3 

_~1, 

+ 13F3I 

,31, 

14  1 

LC3! 

sin 

(* 

+ (0)  - S31 

+ terms  m e 


,ql  '|q  - i(k  + q)e2|  (y*  “ qu) 


(11) 


■Since  e £=  0.015  for  Proton  4,  the  first  of  the  three  terms  within  the  square 
brackets  is  likely  to  be  negligible. 
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3.2  Results 

At  dates  between -14  December  1963  and  6 March  1969,  when  the  orbit  of 

Proton  4 was  appreciably  perturbed  by  the  effects  of  31 ;2  resonance,  there  are 

nine  sets  of  PROP  orbital  elements  available  (orbits  2-10  of  Table  2),  and  ten 

sets  of  US  Navy  elements.  The  19  values  of  inclination  from  these  orbits  are 

plotted  in  Fig. 2.  The  values  were  then  cleared  of  (a)  lunisolar  and  zonal- 

harmonic  perturbations  (with  a combined  maximum  value  of  0.0023°)  and  (b)  the 

effects  of  an  atmosphere  rotating  at  l.Orev/day  (maximum  value  0.0028  ).  The 

PROP  values  of  inclination  were  also  cleared  of  tesseral  harmonic  perturbations 

due  to  the  ^ term  i-n  t^e  geopotential  (maximum  value  0.0018  ).  The  19 

values  of  inclination  after  removal  of  perturbations  are  plotted  in  Fig. 3.  The 

corresponding  observational  values  of  eccentricity  are  plotted  in  Fig. 4,  and 

then  replotted  (on  a much  larger  scale)  in  Fig. 5 after  removal  of  zonal- 

* 

harmonic  and  lunisolar  perturbations. 

The  19  modified  values  of  inclination  and  eccentricity  in  Figs. 3 and  5 
have  been  fitted  with  least-squares  theoretical  curves,  first  separately  and 
then  simultaneously,  using  the  THROE  ’ and  SIMRES  computer  programs  respecti- 
vely, to  obtain  values  of  lumped  31st-order  harmonics  in  the  geopotential. 

First,  the  values  of  i in  Fig. 3 were  fitted  with  equation  (5),  in 
integrated  form,  using  THROE.  The  fitting  was  good,  with  the  measure  of  fit, 
e , having  a value  of  0.96.  This  confirms  that  the  variation  in  i is  consist 
ent  with  a 31:2  resonance  effect,  but  the  numerical  values  of  the  C and  S 
coefficients  are  indeterminate.  They  are  as  follows: 


6-°>2 
io  s31  = 

00 

• 

+1 

• 

o 

,o,: 2 

1°  C3i  - 

0.03 

± 0.8 

6-1’* 

1°  C31  = 

- 87  ± 46 

6-1’1 

10  S31  ' 

- 20 

± 18  | 

► 

(12) 

fi-"1*3 
,0C31  = 

+ 8 ± 12 

fi-"1’3 
10  S31  " 

7 ± 

15 

The  first  two  values  here,  representing  the  (y,q)  = O>0)  terms  in  equation  (5), 
are  (presumably  by  chance)  very  small,  and,  on  substituting  them  into  (5),  we 
find  that  they  comprise  less  than  5%  of  the  total.  These  (y,q)  = O>0)  terms 
are  also  very  likely  to  be  negligible  in  equation  (11)  for  de/dt  ; so  they 
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were  discarded,  and  the  values  of  i were  fitted  with  (y,q)  = (1,1)  and 
(1,-1)  only.  The  resulting  value  of  £ was  0.93,  and  the  values  of  the 
coefficients  are: 


6_ 1 > * 

10°C3]  = - 86  ± 22 


10  C = 8 ± 6 


10  S = - 18  ± 13 


, -1,3 

10°S3]  =6+6 


The  values  are  not  significantly  altered,  but  the  standard  deviations  are 
approximately  halved. 

In  fitting  i , the  sd  of  the  US  Navy  values  was  taken  as  0.003°;  the  sd 
of  the  PROP  values  was  taken  from  Table  2,  except  that  the  sd  for  orbit  2 was 
increased  from  0.0002°  to  0.0005°,  because  of  the  neglect  of  Earth-tide 
perturbations;  and  the  density  scale  height  was  taken  as  50km,  appropriate  to  a 
height  of  280km  (?H  above  perigee).  The  fitted  curve,  shown  as  a broken  line 
in  Fig. 3,  satisfactorily  follows  the  variations  in  the  observational  values. 

Next,  the  values  of  e in  Fig.;<  were  fitted  with  equation  (11),  in 
integrated  form,  using  THROE.  The  (y,q)  = (1,0)  term  was  dropped  - the  values 
obtained  when  it  was  included  being  absurdly  large  and  indeterminate,  as 
expected.  The  measure  of  fit,  e , was  1.48  and  the  values  of  the  lumped 
coefficients  are: 


10  C = - 20  t 15 


10  C3J  = 12  t 4 


10VS31  = 50+11 


c_-  1 > 3 

10°S3]  = 11+4 


In  fitting  the  eccentricity,  the  sd  of  the  US  Navy  values  was  taken  as  0.0001 
(taking  the  sd  as  0.00004  gave  unacceptably  large  values  of  e , though  the 
standard  deviations  of  the  lumped  coefficients  were  not  much  altered);  the 
density  scale  height  H was  taken  as  55km,  appropriate  to  a height  of  325km 
(3H/2  above  perigee).  In  the  course  of  the  analysis,  the  sd  of  e was  doubled 
on  four  of  the  PROP  orbits  (orbits  2,  4,  5 and  10),  since  this  change  improved 
the  fitting.  The  fitted  curve  is  shown  as  a broken  line  in  Fig. 5,  and  it  is 
clear  that  the  US  Navy  values  display  an  oscillation,  of  unknown  origin. 


12 


In  the  US  Navy  elements  as  received,  the  odd-harmonic  perturbations  in  e have 
been  largely  removed  and  have  to  be  restored  (using  the  SDCELS  subroutine).  In 
an  attempt  to  reduce  the  oscillation,  the  restoration  was  made  using  a non- 
standard value  of  , -2.40  x 10  which  includes  an  allowance  for  . 

This  change  reduced  the  amplitude  of  the  oscillation,  but  it  remains  quite 
large. 


The  second  pair  of  (C,S)  coefficients  in  (14)  agrees  well  wiuh  the  second 
pair  in  (13),  but  the  only  agreement  between  the  [l,l]  pairs  of  (C,S) 
coefficients  in  (13)  and  (14)  is  that  both  are  larger  than  the  [-1,3] 
coefficients.  The  discrepancy  between  the  values  of  the  [l,l]  coefficients  in 
equations  (13)  and  (14)  can  be  resolved  by  a simultaneous  fitting  of  the  varia- 
tions in  i and  e . This  has  been  done,  using  the  SIMRES  computer  program^, 
where  there  is  a choice  of  weighting.  Two  alternatives  were  tried:  (a)  with  i 

and  e having  equal  weights,  and  (b)  with  e degraded  by  a factor  of  1.6  in 

recognition  of  its  larger  value  of  e (1.48  as  against  0.93  for  i)  . The 

second  alternative  was  chosen  because  it  is  more  logical  and  because  it  gave 

values  more  consistent  with  (13)  and  (14).  (The  equal  weighting  gave  slightly 
lower  standard  deviations,  however.)  The  fitting  (b),  for  which  e = 1.14  , is 
shown  by  the  unbroken  lines  in  Figs. 3 and  5,  and  the  values  of  the  coefficients 
given  by  SIMRES  are: 


106C3*  = - 45  ± 14 

10  C ■ 10  ± 3 


106S 


1,1 

31 


10°S 


-1,3 

31 


11  ± 7 


2 ± 3 


(15) 


The  values  in  equations  (15)  are  quite  close  to  the  means  of  the  values  in 
equations  (13)  and  (14). 

3.3  Discussion 

The  values  of  the  lumped  coefficients  in  equations  (15)  are  fairly 
consistent  with  those  obtained  from  analysis  of  i alone,  equations  (13),  and 
with  those  derived  from  e alone,  equations  (14)  - they  are  within  about  twice 
the  sum  of  the  sd,  the  {-1,3]  coefficients  being  much  the  closer.  Inevitably, 
the  .simultaneous  fitting  is  a compromise  between  i and  e , and  neither  is 
fitted  so  well  by  the  combined  solution:  the  values  of  s are  1.18  for  i and 
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2.16  for  e , as  compared  with  0.93  for  i and  1.48  for  e in  the  separate 
solutions.  Nevertheless,  the  fitting  of  the  unbroken  curves  in  Figs. 3 and  5 is 
quite  satisfactory,  so  the  combined  solution  is  obviously  preferable. 

The  expected  order  of  magnitude  of  the  lumped  coefficients  can  be 

estimated  by  assuming  that  the  individual  coefficients  of  degree  H have 
5-5  2 

numerical  values  of  order  10  / Z , so  that  those  with  31  ^ l < 43  should 
be  of  order  10  . For  Proton  4 the  numerical  versions  of  equations  (8)  and  (9) 


' C31,3>  - 15C33,3,  + 80C35,3.  ' 2MC37,3I  + 466C39,31  " 534C41,31 


+ 256C,,  + 190C, _ - 328C . n 01  + ...  (16) 

43,31  45,31  47,31 

_-l,3  _ _ _ 

C31  = C31 ,31  " 1 2C33, 31  + 52C35,31  ~ 122C37,31  + 164°39,31  ” 96°41 ,31 


~ 47C43,31  + 104C45, 31  " 8C47,31  + 


and  similarly  for  S , on  replacing  C by  S throughout. 

Since  the  (y,q)  = (1,0)  terms  have  been  dropped,  both  pairs  of  coefficients 
give  lumped  values  of  geopocential  coefficients  of  order  31  and  odd  degree 
(£  = 31,  33,  35  ...).  Equations  (16)  and  (17)  show  that  the  largest  contribu- 
tions are  most  likely  to  come  from  the  coefficients  of  degree  37,  39  and  41. 

Equations  (16)  and  (17)  indicate  that  the  C and  S coefficients  with 
affixes  [l,lj  might  be  expected  to  be  about  3 or  4 times  greater  than  those  with 
affixes  [-1,33.  This  is  in  conformity  with  the  actual  results,  since  ' 
h/ 452  + ll2/\il02  + 22  = 4.5  . 


Equations  (16)  and  (17)  also  suggest  that  if  the  individual  coefficients 
-8  . . . -1,1 
are  of  order  10  , the  expected  magnitude  of  the  lumped  coefficients  Cot  and 

_ — 1,3  _,5  _6  -51 

would  be  about  10  * 10  and  3 x 10  respectively:  the  actual  values 

in  equations  (15)  are  about  three  times  larger.  This  discrepancy  cannot  be 

regarded  as  significant,  since  the  assumed  magnitude  of  the  coefficients 
-5  2 

(10  !%  ) may  be  in  error  by  a factor  of  up  to  2,  and  the  terms  in  equations 

(16)  and  (17)  might  happen  to  add  up  to  considerably  more  than  the  root  of  the 
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sums  of  their  squares.  Nevertheless  these  numerical  values  of  31st-order 
coefficients  should  be  treated  with  the  caution  appropriate  in  any  first  deter- 
mination of  physical  constants. 

4 ATMOSPHERIC  ROTATION  RATE 

The  rotation  rate  of  the  upper  atmosphere  was  estimated  using  15  values 
of  inclination,  given  in  Table  2 by  orbits  11-25,  between  March  &nd  July  1969. 
These  values  were  cleared  of  lunisolar  and  geopotential  perturbations  and  the 
modified  values  are  shown  plotted,  together  with  their  standard  deviations,  in 
Fig. 6 (from  MJD  40293  to  40424).  Similarly  modified  US  Navy  values,  shown  by 
crosses,  are  included  for  general  comparison.  Fig. 6 also  shows  the  values 
during  the  resonance  period  and  the  broken  curve  of  Fig. 3,  after  the  restoration 
of  the  atmospheric  rotation  perturbation. 

The  theoretical  change  in  inclination  was  calculated  for  various  values 
of  atmospheric  rotation  rate  (expressed  as  A times  the  Earth's  rotation  rate), 
using  oblate-atmosphere  theory  , with  numerical  integration  at  intervals  of 
about  6 days  (corresponding  to  22.5°  steps  in  argument  of  perigee). 

A theoretical  curve  fitted  to  the  15  values,  as  shown  in  Fig. 6,  gives  a 
mean  value  of  rotation  rate, 

A = 1.10  t 0.05rev/day  , 

equivalent  to  an  average  zonal  wind  of  40  ± 20m/s  west  to  east,  at  a mean 
height  of  260km,  effectively  averaged  over  latitudes  up  to  about  35°  (since 
atmospheric  rotation  has  little  effect  on  an  orbit  of  inclination  51°,  at 
latitudes  of  35-51°). 

Recent  results  from  analysis  of  many  satellite  orbits11  have  indicated  a 
difference  between  evening  and  morning  winds,  so  the  values  were  separated  into 
two  groups,  with  the  division  where  the  local  time  at  perigee  is  12h  (MJD  40375), 
and  fitted  with  two  curves  as  shown  in  Fig. 7 to  give  morning  (4— 1 2h)  and 
afternoon-evening  ( 1 2— 24h)  values  of  A ; the  local  time  at  perigee  is  marked  at 
the  top.  The  values  of  A obtained,  1.3  ± 0.1  for  afternoon-evening  and 
0.9  ± 0.1  for  the  morning,  agree  with  those  obtained  from  other  satellites  *. 

5 CONCLUSIONS 

Our  analysis  shows  that  the  orbit  of  Proton  4 was  substantially  affected 
by  31:2  resonance  with  the  geopotential  between  December  1968  and  March  1969. 
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The  sets  cT  orbital  parameters  obtained  from  observations  were  accurate  enough 
for  the  variations  in  inclination  and  eccentricity  near  resonance  to  be  analysed 
to  determine,  for  the  first  time,  numerical  values  of  lumped  geopotential 
coefficients  of  31st  order  and  odd  degree. 

The  simultaneous  least-squares  fitting  of  the  values  of  inclination  and 
eccentricity  at  19  epochs,  using  the  SIMRES  program,  leads  to  the  following 
values  of  lumped  31st-order  harmonic  coefficients: 

10°C31  = - 45  ± 14 

, -1,3 

10  C3]  = 10  ± 3 

Equations  (16)  and  (17)  give  these  lumped  coefficients  explicitly  in  terms  of 
the  individual  geopotential  coefficients  of  order  31,  C and  S . 

X/  j J 1 X/  j J 1 

Although  the  numerical  values  of  the  lumped  coefficients  obtained  here  are 

. . 12 
not  of  high  accuracy,  this  work,  and  the  recent  analysis  of  29:2  resonance  , 

show  that  it  should  be  feasible  to  evaluate  individual  coefficients  of  order  31 

and  29  - and  possibly  also  27  and  25  - from  analysis  of  a number  of  resonant 

orbits  at  different  inclinations,  if  enough  accurate  observations  can  be 

obtained  near  the  time  of  resonance. 

The  atmospheric  rotation  rate  A has  been  evaluated  from  the  change  in 
inclination  of  Proton  4 between  Marcn  and  July  1969.  The  mean  value  is  found 
to  be  1.10  ± 0.05rev/day  at  260km  height,  corresponding  to  a zonal  wind  of 
40  ± 20m/s  from  west-to-east,  at  a representative  latitude  near  30°.  Dividing 
the  values  into  two  groups  gave  evening  and  morning  values  of  A , as  follows: 

A = 1.3  ± O.lrev/day  for  afternoon-evening  ( 1 2— 24h)  at  270km  height,  equivalent 
to  a west-to-east  wind  of  120  ± 40m/s;  and  A = 0.9  ± O.lrev/day  in  the 
morning  (4-12h),  for  240km  height,  corresponding  to  an  east-to-west  wind  of 
40  ± 40m/s.  These  three  values  of  A,  which  are  biased  towards  latitudes  of 
0-35°,  agree  well  with  the  recent  findings  from  other  satellites^. 
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Table  2 

ORBITAL  PARAMETERS  FOR  PROTON  4,  WITH  STANDARD  DEVIATIONS 
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